Hole spectral functions in lightly doped quantum antiferromagnets 
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We study the hole and magnon spectral functions as a function of hole doping in the two- 
dimensional (2D) t — J and t — t' — t" — J models working within the limits of the spin-wave 
theory, by linearizing the hole-spin-deviation interaction and by adapting the non-crossing approxi- 
mation (NCA). We find that the staggered magnetization decreases rather rapidly with doping and 
it goes to zero at a few percent of hole concentration in both t — J and the t — t' — t" — J model. 
We find that with doping, the residue of the quasiparticle peak at G = (±7r/2, ±7r/2) decreases 
rapidly with doping and the spectral function is in agreement with high resolution angle-resolved 
photo-emission spectroscopy (ARPES) studies of the copper-oxide superconductors. The observed 
large shift of the chemical potential inside the Mott gap is found to be a result of broadening of the 
quasiparticle peak. We find pockets centered at G, similar to those observed by quantum oscillation 
measurements, (i) with an elliptical shape with large eccentricity along the anti-nodal direction in 
the case of the t — J model, and (ii) with an almost circular shape in the case of the t — t' — t" — J 
model. We also find that the spectral intensity distribution in the doped antiferromagnet has a 
waterfall-like patten along the nodal direction of the Brillouin zone, a feature that is also seen in 
ARPES measurements. 

PACS numbers: 71.10.-w,71.10.Fd,71.27.-f a,74.72.-h,79.60.-i 



I. INTRODUCTION 



During the last two decades there is significant progress 
in angle-resolved photo-emission spectroscopy (ARPES) 
due to improved energy and angle resolution by at least 
an order of magnitudei^— . Furthermore, advances in im- 
proving sample quality and characterization have enabled 
ARPES to become a leading tool in the investigation of 
otherwise hidden information about copper-oxide based 
superconducting (COSC) materials. ARPES directly ex- 
poses the momentum and frequency dependent imagi- 
nary part of the single particle Green's function, a quan- 
tity which sits at the heart of quantum many-body the- 
ory. 

As the wave-number of the incident photon is varied, a 
smooth transition of the intensity resembling a waterfall- 
like pattern is observed in ARPES of the COSC materi- 
als at various hole or electron doping concentrations^"— i^. 
While it has been argued that this could be a matrix 
element effect, similar intensity features are obtained 
by means of theoretical models and recently in differ- 
ent materials and under various conditions^. Quantum 
Monte Carlo calculations showed such intensity pattern 
in 2D Hubbard modelsii. In addition, it was suggestedi^ 
that these waterfall-like features of the ARPES inten- 
sity are caused by the gradual transfer of the spectral 
weight from the lowest energy quasiparticle peak at the 
wavevector G = (±7r/2, ±7r/2) to the higher energy so- 
called "string" states as the incident wave-vector gradu- 
ally moves towards (0, 0) along the nodal direction of the 
Brillouin zone. This suggestion was based on a calcula- 
tion of the single hole spectral function in the t — J model 
by broadening the peaks of the spectral function corre- 
sponding to the string excitations!^. Later, the broad- 



ening required to achieve a reasonable agreement with 
the APRES spectra was obtained by considering the ef- 
fects of optical phonons at finite temperatures^^. In the 
present paper we use the t — J and t — t' — t" — J models 
within the limits of the spin-wave theory by linearizing 
the hole-spin-deviation interaction and within the non- 
crossing approximation (NCA) to study the effects of 
doping on the single-hole spectral functions. We show 
that the broadening of the peaks corresponding to the 
high energy string excitations can be achieved naturally 
as an effect of doping, thus, leading to a smooth waterfall- 
like intensity distribution along the nodal direction of the 
Brillouin zone. 

Sometime ago, ARPES studies^^ on the undoped par- 
ent COSC compounds reveal the existence of a relatively 
well-defined low-energy quasiparticle peak at wavevector 
G which follows a dispersion with a bandwidth of similar 
magnitude to that calculated in Ref. "SG*. Rather recently, 
however, Fournier et aL reported that the quasiparticle 
residue as a function of doping is essentially zero below 
approximately 10% hole concentration. In addition, Shen 
et al.^^ have reported ARPES measurements where the 
quasiparticle peak has a broad line shape (much broader 
than the ARPES resolution) and they also found what 
they call a "paradoxical shift" of the chemical potential 
within the Mott gap. Our results indicate that these 
seemingly contradictory findings can be reconciled using 
the calculations reported in the present paper. We find 
that the lowest energy quasiparticle peak observed near 
(±7r/2, ±7r/2) in the undoped cuprates^'*, is very sensi- 
tive to a small amount of doping. Namely, the quasipar- 
ticle residue rapidly goes to zero as a function of doping 
concentration, and, in addition, it is known that optical 
phonons broaden the peaki^. Furthermore, we find that 
the facts, that the quasiparticle peak is broad and that 
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the chemical potential is shifted relative to the quasi- 
particle peak position, are not two independent effects. 
When we repeated our calculation by broadening the 
quasiparticle peak we found that the broadening itself 
leads to a very significant shift of the chemical potential 
relative to the lowest energy quasiparticle position which 
can be understood in a straightforward manner. 

Furthermore, we find that, at low doping, elliptical hole 
pockets form, which are centered around G, in agreement 
with rather recent quantum oscillation measurements'^®. 
Our calculated intensity of the spectral function shows 
that there are luminous elliptical rings centered around 
the points G of the Brillouin zone. We argue that these 
rings could be related to the banana-like features seen in 
ARPES reported on the under-doped COSC materials^. 
Our formulation, which is based on the existence of 
antiferromagnetic (AF) long-range-order (LRO), cannot 
yield a spectrum sharing the single-banana-like feature 
seen in ARPES, where only one side around the G is 
luminous. The reason lies in the reflection symmetry 
about the line (tt, 0) (0, tt) of the Brillouin zone which 
is present when AF-LRO is present in the system. There- 
fore, we view the banana-like shape as a result produced 
by "scissoring" the elliptical ring along the above sym- 
metry line, an operation which breaks this symmetry re- 
quired by the AF-LRO. The measurements presented in 
Ref. [^were carried out in materials which lack AF-LRO. 

Linear spin-wave (LSW) theory can describe the un- 
doped and lightly doped antiferromagnets as long as AF- 
LRO is present in the system. In the undoped COSC par- 
ent compounds, which arc antiferromagnetic insulators'*^, 
the AF-LRO disappears at a very low doping concen- 
tration. For example, in La2~-xSrxCuOi, the AF-LRO 
is destroyed at doping x ~ 0.015^^. On the other 
hand, magnetic fluctuations persist even at higher dop- 
ing concentrations and short-range spin correlations are 
observed in the inelastic neutron scattering experiments 
at reasonably high doping concentrations^^"— . Calcula- 
tions out — J and related models also show the staggered 
magnetization of the system to go to zero as the undoped 
system is gradually doped with holes^"— . 

In order to study the 2D hole-doped antiferromagnet, 
we have considered the t — J and t — t' — t" ~ J mod- 
els working within the limits of the spin-wave theory, 
by linearizing the hole-spin-deviation interaction and by 
adapting the NCA. Within the NCA the self-energy di- 
agrams are considered for both holes and magnons and 
the Dyson's equation is solved to obtain the one-particle 
Green's function. The contribution of the magnon po- 
larization bubbles are considered and, thus, the magnon- 
mediated hole- hole interaction is included. We should 
mention here that we do not use the rigid band approxi- 
mation which uses the hole bands of the single hole cal- 
culation to obtain the hole spectra at finite concentra- 
tion neither do we use any approximate form of the hole 
spectral function to write it as a sum of a coherent quasi- 
particle peak and an incoherent continuumi^iii. 

Hole and magnon Green's functions for a 2D t — J 



model at finite-doping, using LSW approximation and 
NCA, have been studied earlier (Ref. i27|). However, the 
goal and focus of the present study is quite different; as 
already discussed, there are rather recent high resolution 
ARPES studies which exhibit the existence of the wa- 
terfall features^"— and other studies which reveal the 
shape of the most intense part of the Brillouin zone near 
the Fermi surface-'^. Our calculations are done with the 
aim to understand, if possible, these and related features 
of the ARPES measurements. In addition, as we will 
discuss in the paper, our results differ significantly from 
the results of Ref. [27|. For example, we find that only a 
few percent hole concentration is enough to destroy the 
long-range antiferromagnetic order, while in Ref. [23 re- 
sults for the spectral function were presented for up to 
hole doping x ~ 0.25, where we find that this approach 
breaks down precisely because there is no LRO for this 
level of doping. 

Our approach discussed above relies on the existence 
of AF-LRO both within the model and in the materi- 
als which we plan to apply our findings, therefore, we 
only hope to capture the physics of the real materials 
in the limit of light doping. It is known, however, that 
even in the above models the AF-LRO may be unstable 
upon doping leading to other forms of order, such as spi- 
ral orderiiSi^^i^, stripe ordering^Si^ or phase separation 
instability'^^'™. Therefore, we begin by assuming that 
the AF-LRO is present in the system at low doping and 
we proceed with a perturbative approach which asserts 
that spin wave excitations are well-defined in the system. 
If the AF order ceases to exist and it gets replaced by any 
of the above phases, this should show up as an instability 
in our calculation. We will discuss this further in Sec. lIVI 

The paper is organized as follows: In the next section 
we discuss the formulation and the approximation scheme 
followed in this work. In Sec.|X]of the Appendix we out- 
line some of the technical details of the calculation. In 
Sec. mil we present our results on the hole spectral func- 
tions, the hole energy momentum dispersions and the 
intensity of the spectral functions is compared with the 
waterfall and banana-like features observed in ARPES. 
In Sec. IIV Al our results on the magnon spectral func- 
tions, the magnon dispersion relations, and the staggered 
magnetization as a function of doping are presented. In 
Sec.|V]we summarize our findings and conclusions. 



II. FORMULATION 

The 2D t — J or the t—t' — t" — J model, when expressed 
in terms of operators which create spin-deviations above 
the Neel order and by keeping up to quadratic terms in 
these operators, can be approximated by the following 
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Hamiltoniaii2£ 



The NCA self-energies are given by 



H = 



k k 

^k/k-q[ff(k, q)aq + g{k - q, -q)/3lq] 

k,q 



where /ik, /k are the spin- less hole annihilation operators 
(the so-called slave fermions) for J. and t sub-lattices re- 
spectively and ak,/3k are the magnon annihilation oper- 
ators for the corresponding sub-lattices. Ck, i^k are bare 
hole and magnon energies respectively with efe = for t~ 
J model and = 4t'cos fc^^cos ky + 2t"[cos 2fea;-|-cos 2ky] 
for t — t' — t" — J model. We use the parameters t' = 
-0.33t and t" = 0.22t that gives good fit with ARPES 
results of the energy band along anti-nodal directions'^. 
The coupling constants ^(k, q) are defined in Ref. IH. 




FIG. 1: NCA self-energy diagrams for holes (a,b) and 
magnons (c). 

For a system doped with holes, the Dyson's equations 
for hole / spin- wave excitations become 



G(k,w) = 
D{k,u;) = 



1 



w — ek — S(k, io) + i sgn{uj — fi)r] ' 
1 



uj — fik — n(k, uj) + irj^ 



(2) 



where G'(k, cj) and D(k,uj) are the hole and magnon 
Green's functions and E(k, w) and n(k, w) are the re- 
spective self-energies, fi is the chemical potential for holes 
and rj is the broadening parameter. Here we would like to 
note that for a doping concentration of x in a system of 
total number of sites N, the single-hole Green's function 
describes the spectral properties of a hole (or rather a 
fermionic excitation) added or removed from the system 
where xN number of holes are already present. 

In this paper we will work within the so-called non- 
crossing approximation (NCA), as the higher order ver- 
tex corrections have been shown earlier to be small^i^. 



^(1^'^) = ^E/ rf^[5'(k,q)G(k-q,c.-i.)I?(q,i.) 



/k/^k-q[5(k-q,-q)a'_q + .g(k,q)^q] -I- i/.c. (1) n(k,w) 



V(k - q, -q)G(k -q,uj- v)D{-q, ^v)], 
i 

'2^ 



q 

(3) 



The two terms of the hole self-energy describe the hole 
motion outside the Fermi-sea (of holes) and electron mo- 
tion inside the Fermi-sea respectively. Fig [TJa) describes 
hole motion forward in time where as Fig [IJb) describes 
hole motion backward in time which is equivalent to elec- 
tron motion forward in time. It can be easily shown that 
this second term becomes zero for an undoped system. 

Fig E^c) shows the magnon self-energy or spin- 
polarization bubble. This term is also non-zero only for 
the system with finite doping. 

We can get rid of the infinite limits of the energy in- 
tegral over if we utilize the spectral representations of 
the Green's functions: 



/oo 
de[ 
-OO 



Aik,e) 



+ 



B(k,e) 



Ld — e — ^ + IT] UJ — e — fX — IT] 

C(k,e)de 



^4) 



-oo ^ - e + «'7 



(5) 



where /i is Fermi energy and A(k,e), i?(k,e) are hole 
spectral functions defined as 

S(fc,a;) = ^| (m|/ifclV) P <5(c. + e^-^), (6) 



where {\n)} and {|m)} denote the complete set of ba- 
sis states for [N + 1) and [N — 1) hole/particle systems 
and e^"*"^ and e^~^ are the corresponding energies. Here 
is the interacting ground state of the Hamiltonian. 
G(k, e) is the magnon spectral function defined as 



(7) 



where w„ is the energy of the state \n) which can be 
connected to the ground state via the magnon creation 
operator. 

Using these auxiliary functions, the imaginary part of 
the self-energy takes the following form 
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UJ — fl 



S/(k,w) = dL'g^{k,<i)Giik- q,u} - iy)Diiq,iy)9{uj - 1^ - fi) + / diyg^{k - q, -q) 

1 f 

G'7(k-q,w-i^)7:>/(-q, -i^)6'(/i-Lj + j^)] , n/(fc,a;) = -[/ di^ ^ 5^(q, k)G/(q, w + i/)G'/(q - k, z/) 



The derivation of Eq. [5] involves a convolution of spec- 
tral functions utilizing the following representation for 
the delta function: 6(x) = ^lirrij-i^f) ^2 ^,^2 ■ 

The real parts are obtained by the Kramers-Kronig 
relation 



Sfl(k,a;) 
ni?(k,c^) 



1 



-P 



dcj' 



sgn{uj' — fi)T,i{k,u}'] 
sgn{uj')Ui{k,uj') 



-E 



- 1 



N 

k 



Now, the relationship 

(afc4} = ==z i^(fc,t ^0+) (13) 

(14) 



leads to 



(9) 



(ofca^) = / C{k,e)di 



(15) 



Here, we should mention that though the calculation 
of the imaginary part of the self-energy avoids the infi- 
nite energy integration limits, the calculation of the real 
part by Kramers-Kronig relation does involve an energy 
integration between (—00, 00) (which is done numerically 
in the range {—8t,8t) as discussed in Sec. El- The sum 
over q is avoided in the calculation of the real part. 

Once the Green's functions are obtained we can com- 
pute the sub-lattice/staggered magnetization which is a 
measure of the long range antiferromagnetic magnetic or- 
der in the system. It is given by M5 — S — e where S is 

the spin of the system and ^ — J2k (('^k'^k) + (^k&k)^ , 

the spin deviation and aj^ and fej^. the spin-deviation op- 
erators. For an undoped isotropic antiferromagnet this 
leads to 



(10) 



where 7k = {coskx + cosky)/2. For a 2D spin-1/2 
isotropic antiferromagnet its value is e ~ 0.197. We 
should mention here that the validity of the linear spin 
wave theory is based on the smallness of the parameter 

In the case of the doped system we also have Ms — 
S — e where we find that 

' = ^E[^k((«feaI) + (/3-fe/3^fe))-(?7k + l) 



(11) 
(12) 



i.e., this is the total spectral weight of the magnon spec- 
tral function. 

We are not considering the off-diagonal magnons, i.e., 
(/3_fcQ!fc) — = 0). Thus, the final expression is 

. = l^[-^(2/c(fc,.)d.-l)-l]. (16) 

The technical part of the calculation is discussed in 
Sec. E] of appendix. 



III. HOLE SPECTRAL FUNCTIONS 

A. Spectral functions 

The hole spectral functions for various doping concen- 
trations are shown in Fig. [2] and Fig. [3] (in all of our 
figures where we present spectral functions we do that 
for Lu greater than the Fermi level) for the cases of the 
t — J and t ~ t' ~ t" — J model respectively. The various 
lines in a given graph correspond to different wavevectors 
which are varied along the diagonal (0, 0) (tt, tt) direc- 
tion of the Brillouin zone with the lowest curve corre- 
sponding to k = (0, 0) and the middle to k = (7r/2, 7r/2). 
There is a well-defined quasiparticle peak at (7r/2,7r/2) 
which we have truncated at some finite value in order to 
make the other features of the spectra visible. We can 
see that with increasing doping the quasi-particle peak at 
(±7r/2, ±7r/2) gets truncated at the Fermi level (as it falls 
below the Fermi level upon low doping) and the higher 
energy string states gradually lose their strength and be- 
come broader. Fig. [2] also illustrates that for k = (0,0) 
the lowest energy quasiparticle fades away with doping 
and the higher energy "string" states accumulate most 
of the spectral weight as the doping increases. Namely, 
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FIG. 2: Hole spectral function of the 2D t — J model at 
J = OAt{L = 32) for various doping concentrations x. The 
wavevector is varied along the diagonal (0, 0) —5- (n, n) direc- 
tion of the Brillouin zone. 



starting from the quasiparticle peak at (7r/2,7r/2) and 
moving towards (0,0), we see a gradual transfer of the 
spectral weight from the lowest energy peak to the higher 
energy "string" states. This is exactly what was pro- 
posed in Ref. 12 as the origin of the waterfall- like features. 
Fig. [3] illustrates that a very similar behavior holds for 
the t-t' -t" ~ J model. 



x=0.027 x=0.053 x=0.086 




FIG. 4: Hole spectral function of the 2T) t~t' -t" - J model 
at J = QAt for various doping concentrations x calculated 
with a large value of rj = QAt. 




doping concentration x 



FIG. 3: Hole spectral function of the 2D f - 1' - 1" - J model 
at J = QAt{L = 32) for various doping concentrations x. 
The wavevector is varied along the diagonal (0, 0) (tt, tt) 
direction of the Brillouin zone. 



FIG. 5: Quasiparticle residue for k = (7r/2, 7r/2) as a function 
of hole concentration for the t — J and the t — t' — t" — J model 
at J = OAt. 

Sometime ago Wells et al.— found that there is a well- 
defined quasiparticle peak at (7r/2, 7r/2), within their low 
resolution ARPES data, which follows a dispersion with 
a bandwidth which is also in agreement with the results 
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for the t — J and related modelsiSiSS. Later, however, im- 
proved resohition ARPES raeasurements by Shen et al 
reveal that this low energy peak is significantly broader 
than the limit imposed by the experimental resohition. In 
addition, Shen et al.^^ find that the chemical potential, 
at low doping concentrations x, is significantly shifted 
from the lowest energy quasiparticle peak, thus, appear- 
ing to enter the Mott gap. Furthermore, rather recently 
Fournier et ali^ find that the quasiparticle residue is es- 
sentially zero below doping concentration of 0.1. Next 
we will try to understand these results within the re- 
sults of our calculation. First of all, our present calcu- 
lation which ignores the role of phonons shows that the 
lowest energy quasiparticle peak at G is much sharper 
than the experimental peak. Previouslyii, we added the 
role of optical phonons and of finite temperature, both 
of which have been neglected in the present calculation, 
and we demonstrated that the lowest energy peak be- 
comes broadened. In Fig. |4] we present the results of our 
calculation obtained using the t — t' — t" — J model for 
three different values of doping x and a relatively large 
value of the broadening parameter 77, in order to com- 
pare with Fig. 3 of Ref. [l^ Notice that the evolution of 
spectra seem qualitatively very similar to those reported 
in Ref. ITsl . In particular, we notice that the broadening 
causes the position of the chemical potential (the value of 
u! in the spectra in Fig.|4]has been shifted by the value of 
fi such that the position of /i is at = 0) to move signif- 
icantly away from the quasiparticle peak inside the Mott 
gap. This is very similar to the observationiS,. Therefore, 
we find that the results reported in Ref. [l^ regarding 
what was referred to as "paradoxical shift" of the chemi- 
cal potential, can be understood as a consequence of the 
broadening. 

Fig.[5]illustrates the quasiparticle residue as a function 
of doping for the t — J and the t — t' — t" — J model for 
J/t = 0.4 obtained by integrating the spectral function 
from w = /i up to the first minimum of the spectral func- 
tion for uj > ji. Our results on the t — t'~t" — J model in- 
dicate that the quasiparticle residue becomes very small 
away from the undoped limit and as can be inferred from 
Fig. 4 of Ref. |^ our results are within error bars of the 
experimental values. In addition, Fournier et al.^ provide 
no data in the immediate vicinity of zero doping. 

Therefore, we can conclude that our results can explain 
the ARFES^iiiii^ results without additional hypotheses. 



B. Hole energy dispersions 

Fig. |6] shows a 2D contour plot of the rigid hole energy 
dispersion E{kx,ky) derived from the t — J model on a 
32 X 32 lattice for J — OAt. As we follow the iso-energy 
lines we can see that the low energy contours are ellipses 
around the wave vectors (±7r/2, ±7r/2) (the lowest en- 
ergies are indeed at (±7r/2, ±7r/2) which is not shown 
in the plot). At higher energy, the topology gradually 
changes and the ellipses increase in size and become con- 




FIG. 6: Contour plots of the rigid hole energy dispersion 
E{kx, ky) in the t ~ J models for J = OAt in a 32 x 32 lattice 



nected to each other giving rise to new pockets around 
(0,±7r) and (±7r, 0). This transition occurs at doping 
around x ~0.145. At even higher energy these newly 
formed pockets disappear and disconnected contours are 
obtained centered around (0,0) and (±7r, ±7r) and hap- 
pens above x ~ 0.19. A similar change of topology as 
function of x occurs in the t — t' — t" — J model; the main 
difference being that the fc-space anisotropy of the spec- 
tral intensity is reduced making the hole pockets more 
circular than elongated. 

The contour plot of the energy dispersion obtained for 
doped system show similar features at low doping. The 
hole energy dispersion E{kx,ky) is obtained as the low- 
est energy peak above /i in the spectral function as a 
function of k. Fig. [7] shows the contour plots for hole 
energy dispersions at various dopings. We can see that 
the Fermi surfaces that are small hole pockets around 
(±7r/2, ±7r/2) for very small doped system, gradually get 
elongated along anti-nodal directions as the system is 
doped with more and more holes. 

In Fig. ISl the calculated hole spectral intensity plots 
at Fermi energy, for the t — J model (for x = 0.084, top 
left), and for the t-t' -t" -J (for x = 0.094, top right) 
are compared with the experimental ARPES plot (bot- 
tom) for the Na2-xCaxCu202Cl2 sample with x — 0.1 
taken from Ref. 2. The experimental intensity plot has 
a banana-like shape close to (±7r/2, ±7r/2) at this dop- 
ing. Our numerical results show elliptical rings centered 
about (±7r/2, ±7r/2) stretched along off-diagonal direc- 
tion for the t — J model. Our calculations are restricted 
in the regime where there is long-range antiferromag- 
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-3-2-10123 



FIG. 7: Contour plots of the hole energy dispersion E{kx, ky) 
for the t — J model for J = OAt in a 32 x 32 lattice. Top: 
fi = -2.20t, X = 0.022; middle: fj. = -2.10t, x = 0.082; 
bottom: = -2.02t, x = 0.118. 




n 



FIG. 8: Intensity plots of hole spectral functions at the Fermi 
energy for J = OAt as calculated on a 64 x 64 lattice for 
the t — J model (top left) at doping concentration x — 0.084 
and for the t — t' ~ t" — J model (top right) at doping con- 
centration X = 0.094. Bottom: ARPES intensity plot of 
Na2-xCaxCu202Cl2 at X — 0.1 done in Ref. 3 (this part 
of the figure is taken from Ref. '3^ which was reproduced from 
Ref. 3 and the red elliptical Fermi arc was added). All plots 
show only the 1st quadrant of the Brillouin zone. 




(0,0) (71,0) (ii,it) (i;2,ii;2) (0,0) 



FIG. 9: Hole energy band in the t — J model for J — OAt in a 
32 X 32 lattice along (0,0) (tt, 0) {tt,-k) (7r/2,7r/2) 
(0, 0) cut in fc-space. 
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netic order present in the system where the hole spec- 
tra should be symmetric against reflections with respect 
to the (0, tt) (tt, 0) line in the Brillouin zone. As a 
result spectra having banana-like shape cannot exist in 
the t — J or t — t' — t" — J model when there is AF or- 
der. It is possible when the long-range antiferromagnetic 
order is destroyed upon further doping, where the reflec- 
tion symmetry with respect to the (0,7r) — >■ (7r,0) line 
in the Brillouin zone is broken, that the ring-like shape 
of the spectral intensity plot might transform to half el- 
liptical ring-like shape which might resemble the banana 
shape. For example, if we cut the intensity distribution 
obtained for the t — t' — t" — J model shown in Fig. [8] 
(top right) along the (0, tt) -> (tt, 0) line, one of the pieces 
will resemble the experimental intensity plot. The next 
nearest neighbor hoppings present in the t — t' — t" — J 
model reduce the fc-space anisotropy of the hole pockets. 
However, the conclusion drawn from quantum oscilla- 
tion measurementsiS. indicate that there may be elliptical 
pockets centered around the (±7r/2, ±7r/2) points. These 
are shown by the ellipse in Fig. [8] (bottom) . Therefore, it 
is not clear why only one side of the elliptical Fermi arcs 
are luminous in the ARPES studies. The ARPES result 
indicates that the quasiparticles are much better defined 
only just outside of one side of the Fermi arcs. 

We should mention here that because of the finite en- 
ergy resolution rj in our calculation, in Fig. [5] we have 
also smoothed the spectral functions as follows 

A{k,oj) = [ A{k,oj')— -div' (17) 

with T] = 0.02i. 

Fig. ini shows the hole bands for the t — J model along 
(0,0) -5> (7r,0) -> (7r,7r) ^ {71/2,1: /2) (0,0) direc- 
tion in /c-space. The reason that a considerably larger 
bandwidth of the order of t is found at higher doping 
{x = 0.118) (also shown in Ref. \2% is that the lowest 
energy string state, as k approaches (0, 0) or (tt, tt), loses 
all of its spectral weight which is transferred to higher 
energy string-states. However, the staggered magnetiza- 
tion {Ms) is zero at those doping concentrations as can 
be seen from Fig. [19] and so our calculation of the hole 
band for this value of x is outside the range of validity of 
our calculation. 



C. Waterfall-like features 

The single-hole spectral function as obtained from the 
t — J model has high energy string states'^® and these 
states, along with the quasi-particle peak, broaden in 
the presence of optical phonons and finite temperature^'^. 
This enables a smooth transition of intensity from the 
low energy peak at (7r/2, 7r/2) to the high energy peak at 
(0,0)i^ as seen in the experimental waterfall-like feature^. 
Here, we find a similar fading away of the higher en- 
ergy string states with doping and that makes a bridge 
between the high energy and low energy states so that 




(0,0) kp (0,0) kp (0,0) kp 



FIG. 10: Intensity plots in a 32 x 32 lattice from a 2D t — J 
model for doping concentration x = 0.02 (left) and x = 0.09 
(middle,right) for J = QAt. Lorentzian broadening with rj = 
0.03f is used in plot at right. 




FIG. 11; Waterfall feature: ARPES intensity plot^ of under- 
doped Bi2212 (left) and intensity plots in a 64 x 64 lattice 
from a 2D t — J model for doping concentration x — 0.084 
(middle) and from at — t' — t" — J model for x = 0.094 (right) 
for J = OAt. 
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the continuous change in the intensity gives rise to a 
waterfall-like pattern. 

The spectral functions are smoothed as discussed by 
means of Eq. [17] to smooth out the "sharp edges" pro- 
duced by the combination of the finite-size system and 
the finiteness of the energy resolution used in the numer- 
ical calculations. 

Fig. [TU] shows intensity plots along the nodal direc- 
tion obtained using the t — J model for a; = and 
X — 0.09. The discreteness of the high intensity close 
to Fermi wave- vector fc^^^ is due to finite size effects and 
can be smoothed out when using higher values of 77. Fig. 
111! compares intensity plots obtained from the t — J and 
the t — t' — t" ~ J model for a 64 x 64 size lattice with 
the ARPES waterfall-like intensity plots on the under- 
doped Bi2212. Notice a smooth transfer of intensity in 
the t — J model intensity plots resembling the experimen- 
tal ARPES plot^. 
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FIG. 13: Magnon spectral functions for along the (0, 1) di- 
rection of the Brillouin zone obtained for the t — J model for 
J = OAt (L=32) for x = 0.022. 



IV. MAGNETIC ORDER 



A. Magnon spectral functions 



o^x=aoii 




m/t 



FIG. 12: Magnon spectral functions for along the diagonal of 
the Brillouin zone obtained for the t — J model for J = QAt 
(L=32) for X = 0.022. 

In Fig. [Inland Fig. [T3]we present the magnon spectral 
function for x = 0.022 (for a 32 x 32 size lattice) of the 
t — J model along the diagonal and the (1,0) directions 
of the Brillouin zones respectively. For both the undoped 
and doped system, D{k,uj) has a rather featureless peak 
at the renormalized magnon energy flr{k) with a width 
of similar order of magnitude to flrik) which is of order 
J, unlike the hole spectral function, which has a spread 
from approximately -4t to 5t in energy. The renormalized 
magnon energy decreases rapidly with increased doping 
and at some low value of doping concentration x the spin 
wave spectrum collapses at small wave-vectors and this 
signals that the AF-LRO disappears (see also Ref. Isol ). 



G-O x=0.005 
Q-O 0.030 
C^'' x=0.051 



03 - 




(0,0) (0,ji) (ir/2,n/2) (0,0) 



FIG. 14: Magnon frequency in the t — J (top) and t — t' — t" — J 
(bottom) models for J = 0.2t. 
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Fig. [14] shows the spin wave dispersion as a function of 
doping obtained from the t — J and t — t' — t" — J models 
in a 32 X 32 lattice. Doping lowers the magnon frequency 
and above a certain value of x it becomes zero starting 
at low wave- vectors. 
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FIG. 15: Magnon spectral functions for along the diagonal of 
the Brillouin zone obtained for the t — J model for J = OAt 
(L=64) for X = 0.084. 
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FIG. 16: Magnon spectral functions for along the (0, 1) di- 
rection of the Brillouin zone obtained for the t — J model for 
J = 0.4f (L=64) for x = 0.084. 

Fig. [Island Fig.[T6]shows the magnon spectral function 
for X = 0.084 (for a 64 x 64 size lattice) of the t — J model 
along the diagonal and the (1,0) directions of the Bril- 
louin zones respectively. Notice that in this case however, 
along with the broad magnon peak at the renormalized 
magnon frequency there is a sharp peak at = for 
small wavevectors. This is an indication of an instability 
which was also found by Krier et al4S. This instability 
could be due to the spiral orde r^^i^^i'^^ or due to stripe 
ordering^S. 



B. Staggered magnetization 



0.3 - 




0.02 0.04 0.06 0.08 



FIG. 17: Staggered Magnetization in the t—J model for L=20 
and 32 (J = 0.2t). 




X 



FIG. 18: Staggered Magnetization in the t — J and t—t' — t" — J 
models for J = Q.2t (L=32). 

The staggered magnetization [Ms) for a spin- 1/2 2D 
antiferromagnetic square lattice is 0.30341. with hole 
doping, Ms decreases from the above value. Fig. [17] in- 
dicates that the role of finite-size effects may be insignif- 
icant. We find that Ms decreases steadily with doping 
and eventually goes to zero within a doping range be- 
tween 0.07 and 0.08 for the value of J = 0.2i (as seen 
in Fig. [TTI and [18)) . Using larger values of J increases 
the value of the critical doping {xc) where Ms ceases to 
exist. Fig. [T^ shows that for J = 0.4i the magnetization 
goes to zero at a concentration Xc ^ 0.12. Our stag- 
gered magnetization plots are similar to that obtained 
in Rcf. [2l[ - [25l Following the discussion regarding the 
possible instabilities of our previous section, we would 
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0.02 0.04 0.06 0.08 0.1 0.12 



FIG. 19: Staggered Magnetization in the t — J model for 
J = Q.2t and QAt (L=32). 

like to comment that our calculation cannot definitely 
distinguish the case where the antiferromagnetic order 
disappears due to the presence of the stripe or spiral or 
phase separation instabilities3-- at any doping from the 
case where the order disappears at a finite concentration. 
However, the fact that the staggered magnetization feels 
weak finite-size effects is an indication of the latter sce- 
nario. 



V. CONCLUSIONS 

We have studied the 2D t — J model and the related 
t — t' — t" — J model in the low doping regime where the 
antiferromagnetic LRO is present. Our approach requires 
the existence of AF-LRO in order to be qualitatively cor- 
rect because it starts by transforming these Hamiltonians 
by keeping up to quadratic terms of operators which cre- 
ate spin-deviations above the Neel state. We find that 
doping quickly reduces the degree of magnetic order and 
completely destroys it for doping concentration greater 
than X ~ 0.08 - 0.1 (for 0.2 < J/t < 0.4). Experimental 
findings on cuprate superconductors also show quick dis- 
appearance of long range magnetic order with dopinglS^ 
and when the LRO is destroyed by doping the correlation 
length is several lattice constants. 

We find that the residue of the quasiparticle peak at 
(±7r/2, ±7r/2) dies out quickly with doping and in par- 
ticular in the t — t' — t" — J model it becomes negligi- 
ble away from the undoped limit. In addition, we find 
that when we broaden the lowest energy quasiparticle 
peak, the spectral function shares most of its features 
with those obtained by high resolution ARPES studies^^. 
In particular, we find that what is called in Ref.[ll "para- 
doxical shift" of the chemical potential within the Mott 
gap, is a natural result of the broadening of the lowest 
energy quasiparticle peak. Our results indicate that the 
high resolution ARPES results obtained by Shen et alJ^ 



and those recently reported by Fournier et in the low 
doping regime can be understood in a rather simple way. 

Furthermore, we studied the dependence on hole 
doping of the hole pockets which start forming at 
(±7r/2, ±7r/2) at very low hole doping when antiferro- 
magnetic order is present in the system. We find that 
as we increase doping the Fermi "surface" starts as four 
ellipses around the wave vectors (±7r/2, ±7r/2) and the 
topology gradually changes where the ellipses increase 
in size and become connected to each other giving rise 
to new pockets around (0, ±7r) and (±7r,0). At even 
somewhat higher doping these newly formed pockets dis- 
appear and disconnected parts of the Fermi surface are 
obtained centered around (0,0) and (±7r, ±7r). A sim- 
ilar change of topology as function of x occurs in the 
t — t' — t" — J model, the main difference being that the 
fc-space anisotropy of the spectral intensity is reduced 
making the hole pockets more circular than elongated. 
Similar pockets have been inferred from quantum oscil- 
lation measurements^ at relatively low doping. 

We also find that the intensity plot of the spectral 
intensity has features with a ring-like shape which are 
centered around the wave vectors (±7r/2, ±7r/2). Our 
calculations are restricted in the regime where there is 
long-range antiferromagnetic order present in the system 
where the hole spectra should be symmetric against re- 
fiections with respect to the (0, tt) (tt, 0) line in the 
Brillouin zone. As a result spectra having banana-like 
shape as seen experimentally^, cannot exist within the 
models studied here when there is AF order. We argued, 
however, that it is possible when the long-range antifer- 
romagnetic order is destroyed upon further doping that 
these features with the ring-like shape of the spectral in- 
tensity plot could transform to a half ring-like shape (i.e., 
a ring centered at (7r/2,7r/2) which is scissored along the 
(0, tt) — i> (tt, 0) line) which might resemble the banana 
shape pockets near (7r/2,7r/2) seen in the experiment^. 

We have also examined the role of doping in produc- 
ing the waterfall-like features in the spectral function as 
seen in ARPES. In particular, it was suggested in Ref. [l^ 
that these waterfall-like features is the result of spectral 
weight transfer from the lowest energy quasiparticle peak 
to the higher "string" states as the momentum varies 
from (7r/2,7r/2) to (0,0) along the diagonal of the Bril- 
louin zone. We find that doping significantly broadens 
the high energy "string" states and this allows a gradual 
and continuous fiow of spectral weight from the quasi- 
particle peak at (^, ^) to the high energy peak at (0,0). 
We also find that doping does not broaden the lowest en- 
ergy peak enough to have agreement with the experimen- 
tal waterfall-like feature, thus, other hole decay mecha- 
nisms such as decay into phonons or broadening due to 
finite temperatures may be needed. 
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Appendix A: Technical part 

For the numerical calculation, we always keep the en- 
ergy resolution, i.e., the finite-frequency step Aw (O.OOSt) 
less than rj (we used 77 between 0.4t and O.Olt) so that the 
Lorentzian feature of the bare Green's functions can be 
represented properly as also discussed in Refs. Il3ll36l . 
The energy range between which we iterate both the 
Green's functions are {—8t, 8t) as the variance of the spec- 
tra are always well within such limit. We should mention 
here that the Auj used in Ref. [l^ that gives 400 equally 
spaced points in the range -5t to 4t is larger than the 
Lorentzian broadening, rj = O.OlSt that they used. 

We find that the convergence of the self-consistent 
equations is slow for the doped antiferromagnet (com- 
pared to the calculations done in Ref. 36 for obtaining 
the hole Green's functions from undopcd system). Com- 
pared to the t — t' — t" — J model the convergence is even 
slower for the t — J model as they have a complete degen- 
eracy in the bare hole energy. In order to obtain a quick 
convergence we update the spectral functions, in each it- 
eration step, by mixing the values obtained from the last 
two iterations- a technique also adopted in Ref. H^- We 
also notice poor convergence of self-consistent equations 
as Ms becomes very small. 

When calculating the staggered magnetization, we use 
Ci{k,uj) — —D^j'^\k,Lu)/'K for both positive and nega- 
tive ui. Numerically a Lorentzian of width rj is used to 
describe the imaginary part of the bare Green's function 
d'^q^\ which actually should be a delta function at the 
non-negative energy J7fc. And, thus, to obtain the nor- 
malization correctly, we also need to consider the tail of 
the Lorentzian even though that extends to negative oj. 

We note that in calculating the imaginary part of 
the hole self-energy, the energy convolution is utilized 
to get rid of the infinite integration limits. However 
in that derivation the Lorentzian approximation of the 
delta function overlooks the finite Lorentzian- width com- 
ing from the finite 77. 

We find that the density of states below the quasi- 
particle peak, as obtained from the numerical calcula- 
tion, depends strongly on the value of the broadening 
parameter, 77 used. This is an error which comes as a 
result of the finite 77 used in numerical calculations but 
this 77-dependence reduces significantly as higher values 
of Fermi energy (i.e., higher doping concentrations) are 
considered. Table H] illustrates the result for J = 0.4i on 
a, t — J model rigid band. The numerical results coming 
from a system where the Fermi level stays much higher 
than the quasi-particle energy has a smaller error due to 
finiteness of the 77-used. While studying the hole energy 
bands from the doped system we have, thus, considered 
only the cases where the Fermi energy is higher than the 
quasi-particle peak position (which is at ~ — 2.21i for 
J = 0.4t). 

The density of states (DOS) at low energies (close to 
Fermi energy of the undoped system) becomes modified 




FIG. 20: Hole density of states (number of states in each of 
our 5uj) in the t — J (top) and t — t' — t" ~ J (bottom) models 
for J = 0.2t (L=32). The positions of jj, are also shown. 



Q—Q Doped t-J model results 
B—B Rigid band l-J model res 




FIG. 21: Variation of fi with doping concentration x for J = 
OAt in a 32 X 32 lattice. The blue star on the y-axis at energy 
-2.21t shows the quasi-particle peak position of the undoped 
antiferromagnet. 
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TABLE I: Doping concentrations at different /i for different sets of rj 





-2.7 


-2.5 


-2.4 


-2.3 


-2.2 


-2.1 


-2.0 


0.1 


0.0028 


0.0041 


0.0053 


0.0080 


0.0269 


0.142 


0.216 


0.01 


0.0019 


0.0027 


0.0034 


0.0053 


0.0218 


0.145 


0.217 



with doping (Fig. [501) the DOS low energy peaks 
move towards higher values at higher dopings (also see 
Ref. [27I ). This shifting of the low energy states towards 
higher energies causes the chemical potential to change 
rapidly with doping as compared to that calculated on a 
single hole rigid band as shown in Fig. [21] In the range 
of doping in which we are able to study the doped system 
within the realm of the LSW approximation we do not 
see fi to cross the DOS low energy peak for the t — J 
model. In the t — t' — t" — J model, on the other hand, 
the low energy peaks in the DOS are much flatter and /i 
is seen to cross a few of them even at low doping (Fig. 

As was discussed in the introduction our results differ 
significantly from those of Ref. [13. For example, we find 
that the staggered magnetization decreases very rapidly 
with doping becoming zero for x ^ 0.08 (at J/t — 0.2). In 
contrast, in Ref. [27l results were presented for much larger 
values of doping concentration x (~ 0.25), where accord- 



ing to our calculation, this type of calculation, which 
requires presence of antiferromagnetic long-range order, 
is invalid. The dependence of our magnon spectra with 
doping is also consistent with the vanishing of the stag- 
gered magnetization at x ~ 0.08. Namely, the renormal- 
ized magnon energy softens with doping, for small k val- 
ues, becoming practically zero for x ^ 0.08 (for J — Q.2t). 

Because of the above disagreement we subjected our 
code to several tests. We calculated the self-energies for 
a single iteration using our program and the result agrees 
with the results which we obtain using the mathematical 
expressions obtained analytically. Second our results for 
the spectra at very small doping matches well with the 
single hole spectra (as computed in Ref. l36h at energies 
higher than the quasi-particle peak of the undoped an- 
tiferromagnet. We also check the dependence of energy 
resolution (Aw) on our numerical integration (Ialj) and 
find that / — /aw ~ (Aw)^, as it should be. 
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